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Abstract. We establish rigorously convergence of a semi-discrete upwind 
scheme for the nonlinear variational wave equation Utt — c(u) (c(u)u x ) x = 
with u\t=o = "0 and ut\t—o = Vo- Introducing Riemann invariants R = 
ut + cu x and S = ut — cu x , the variational wave equation is equivalent to 
R t - cR x = c(R 2 - S 2 ) and St + cS x = -c(R 2 - S 2 ) with c = c' /{Ac). An 
upwind scheme is defined for this system. We assume that the the speed c is 
positive, increasing and both c and its derivative are bounded away from zero 
and that i?|t=o> S\t—o £ L 1 (WL) n L 3 (WL) are nonpositive. The numerical scheme 
is illustrated on several examples. 



1. Introduction 
In this paper we consider the nonlinear variational wave equation 

__ c(u) _^ (M) _j =0 , 
Ou 

u(x,0) = u (x), — (x,0) = v (x), 
at 



(1.1) 



in the strip (x, t) eU T =Rx [0, T). 

The equation, which can be derived as the Euler-Lagrange equation for the 
variational principle 6 ff(ipf — c 2 (tp)ipl) dxdt — 0, can be used to model liquid 
crystals, see Consider namely a nematic liquid crystal in the regime where 

inertial effects dominate. In that case the liquid crystal can be described by the 
director field n = n(x, y, z, t) £ K 3 with ||n|| = 1 that describes the direction of the 
elongated molecules that constitute the liquid crystal. Its potential energy density 
is described by the Oseen-Franck functional 

W(n, Vn) = a \n x (V x n)| 2 + /3(V • n) 2 + 7(n • Vn) 2 , 

where the constants a, (3, and 7 describe the liquid crystal. The dynamics of the 
director field is given by a least action principle 

S f 

— (n t -nt-W(n,Wn))dxdydzdt = 0. (1.2) 
ou J 

Consider the class of planar deformations given by 

n = cos(u(x, t))i + sin(u(x, t))j (1.3) 
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where i and j are unit vectors in the x and y direction, respectively. In that case 



the least action principle (1.2) reduces to (1.1 1 with 



c 2 (it) = a cos 2 u + /3sin u. 



We here analyze ( 1 . 1 1 with more restrictive assumptions on c, as is done in the 
mathematical literature, namely that c is positive, strictly increasing and bounded 
away from zero. We note that ( |1.1[ ) is closely related to the Hunter-Saxton equation, 



which is obtained by a further asymptotic expansion of (1.1), see [B]. 

While short-term existence of regular solutions follows by the Kato method, it 
is clear that the solution in general develops singularities in finite time, even from 
smooth initial data, see HI 131. 



In a series of papers, Zhang and Zheng jTTJ [TU [T3J HU [TS] have analyzed ( 1.1 ) 
carefully, using the method of Young measures. From their many results we quote 
the recent one 14. Thm. 1.1] where they show existence of a global weak solution 
for initial data uo G H 1 (M.) and vq G L 2 (M). The function c is assumed to be 
smooth, bounded, positive with derivative that is non-negative and strictly positive 
on the initial data uq. Their results, and also the relationship to the Hunter-Saxton 
equation are surveyed in |15j . The uniqueness question is open. 

Another approach to the study of ( |1.1| was recently taken by Bressan and Zheng 
[lj. Instead of following the approach based on Young measures, they rewrite the 
equation in new variables where singularities disappear. They show that for uq 
absolutely continuous with Uq x ,Vq E L 2 (R) the Cauchy-problem (1.1) allows a 



global weak solution with the following properties: The solution u is locally Holder 
continuous with exponent |, and the map t h- ► u(t, ■) continuously differentiable 
with values in L^ oc (M) for 1 < p < 2. Further properties are obtained, in particular, 
it is shown that the associated energy, treated as a measure, is conserved in time. 
Up to now, little has been known about the behavior of numerical schemes for the 



equation (1.1 1. Except for some numerical computations in [3J, there are, to the best 
of our knowledge, no rigorous results regarding any numerical method for ( |1.1[ ), and 
the main purpose of this paper is to remedy this situation. Here we introduce a semi- 



discrete upwind scheme for the initial-value problem (1.1 1, i.e., a finite difference 
approximation of the spatial variation, keeping the continuous temporal variation. 
For this scheme we show convergence to a weak solution of ( |1.1[ ), and thus this 
proof offers a different existence proof compared with the others. In addition it 
provides a constructive approach to the initial- value problem in the sense that the 
difference scheme supplies a numerical tool to compute the solution, see Section [4] 
Indeed, we show how the difference scheme performs, both on examples where the 
scheme is proved to converge and otherwise. A similar analysis has been applied to 
the Hunter-Saxton equation, see [5]. 

We now turn to a more detailed and technical discussion. Weak solutions are 
defined as follows. 

Definition 1.1. Let denote the set R x [0, T), T > 0. By a weak solution u of 



pTTj ) we mean a function u E L°°([0, T]; W 1:P (M.)) <~)C(Ht), u t G L°°([0, T];L P (E.)), 
for all pe [1,3 + q], where q is some fixed positive constant q > 0, such that 

J J (pt<p dtu — (?{u)d x ip d x u — d (u)c(u)ip (d x u) 2 ) dxdt = (1-4) 



for all test functions ip E C^°(Ht)- The initial values are taken in the sense that 
?([0,T];L 2 

n-r when £ — * OH 



u( ■ , t) — ► no in C([0, T); L 2 (R)) as t — > 0+, and d t u{ ■ , t) ~ * vq as a distribution in 
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A common approach to ( 1.1 1 is first to re- write the equation in terms of Ricmann 
invariants. To that end we define 



du , s du „ du , s du 

and the auxiliary function 



4c(u) 



(1.5) 



Then the wave equation (1.1) is formally equivalent to the 3x3 system 
R t - c{u)R x = c(u) (R 2 - S 2 ) , 
St + c(u)S x = -c(u) (R 2 - S 2 ) , 

u x = irr ^{R-S), 
2c(u) 

Ro = R\t=o = v + c(u )u' Q , S = S\ t =o =v - c(u )u' . (1.6) 
Clearly we also have 

Ut=*(R + S). (1.7) 
In order to make this well defined, we use the boundary condition 



(1.8) 



lim u(x, t) = 0. 

x — > — oo 

The equations for R, S can also be written on conservative form, viz. 

Rt - {c(u)R) x = -c(u) (R-S) 2 , 

S t +(c(u)S) x = -c(u) (R-S) 2 . 

Throughout the paper we will assume that c is a Lipschitz continuous function such 
that 

< d < c(u) < C 2 , and < c'(u) < M x . (1.9) 

The approach by Zhang and Zheng based on Young measures follows two distinct 
routes. Either one can use a viscous regularization of the system (1.5 1 by adding 
the terms eR xx and eS xx to the first and the second equation, respectively, and 
subsequently analyze in detail the behavior of the solution as e — > 0, see [12] . 
Alternatively [TTJ [T^l E] , one can replace the quadratic growth on the right-hand 
side of equation ( |1.5[ ) by a linear growth for large values of R 2 and S 2 . More 
specifically, introduce the function 



Qe(P) 



'f(|P|-i) for|P|>i, 
P 2 for IPI < i, 



and replace the terms R 2 and S* 2 by Q € (R) and Q e (S), respectively, in the first 
and the second equation. Again the behavior of the solution has to be analyzed 
carefully as e — > 0. 

Our approach is based on Young measures for a semi-discrete finite difference 
upwind scheme. More precisely, introduce a positive discretization parameter Ax, 
and approximate R(jAx,t) and S(jAx,t) by Rj(t) and Sj(t), respectively, that is, 
R(jAx,t) « Rj(t) and S(jAx,t) ~ Sj(t), j £ Z. Observe that we keep the time 
variable continuous. The dynamics of (Rj(t), Sj(t)) is governed by the system of 
ordinary differential equations 

R' (t) - c^ +1/2 (t)D + ^(t) = Cj(t) (R 2 (t) - S 2 (t)) , 

S' 3 (t) + c,_ 1/2 (f)D-S,(t) = -cj(t) (R 2 {t) - S 2 {t)) , 
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where D±Kj = ±(Kj±i — Kj)/Ax. Furthermore, the functions Cj±i/2 and Cj are 
defined as functions of Rj and Sj, cf. Section [2] The system is augmented by 
appropriate initial-data. We recover the function ua x by the formula 

2c(u)du= / (RAx(x,t) - SAx(x,t)) dx, 



o J 
where Rax equals Rj(t) on [(j — |)Ax, (j + |)Ax), and similarly for Saz- We first 
show that the system possesses solutions that are local in time, and a subsequent 
a priori estimate turns the local solution into a global one. Once the existence of 
solutions of the ordinary differential equations has been established, we follow the 
approach of Zhang and Zheng closely. 

Formally, a smooth solution of ( |1.8| will satisfy the identity 

(f(R) + f(S)) t - (c(u)(f(R) - f(S))) x = 2cH(R, S), (1.10) 

where 

H(R, S) =\{R 2 - S 2 )(f(R) - f(S)) - (f(R) - f(S))(R - S), (1.11) 



for any smooth function /. The corresponding discrete relation, Lemma |3.1[ is 
rather more complicated. However, based on this, one shows that the difference 



scheme keeps the L 2 norm of {Rj,Sj} from increasing, cf. Corollary 3.2 a similar 
result holds in the continuous case as well, cf. [T2J Lemma 1]. Intrinsic to the 
equation is a blow-up property that is not fully understood. Indeed it is known, 
see |4j , that there exist examples with R and S of opposite sign initially, where the 
solution becomes unbounded. However, if the initial data both are negative initially, 
one can show that the solution remains regular, see, e.g., |12l Thm. 2]. Henceforth 
we will make the assumption here that R and S are nonpositive initially. As in the 



continuous case, [12 Lemma 4], one can show also in the discrete case, Lemma 3.3 
that the equation enjoys invariant domains: If (Rax, Sax) both are nonpositive 
initially, then they will remain so. Furthermore, if (Raxt Sax) m addition are 
bounded from below initially, they will remain so, with the same lower bound. 
From this it follows that L p norms do not increase, cf. |12l Lemma 5] and Lemma 



3.4 Using this one can show in the discrete case, cf. Lemma [376] using the Arzela- 
Ascoli theorem, that there exists a function u such that 

uax - > u uniformly on compacts in K x [0,T]. 

The remaining part of the analysis is to show that the limit indeed satisfies the 
equation. From a priori L p bounds we infer that Rax R and Sax ^ S in 
L°°([0,T];L 2 (R)) (recall that R Ax equals Rj(t) on [(j - ±)Ax,(j + ±)Ax), and 

similarly for Sax), and ( Rax - Sax) 2 (Rax - Sax) 2 in L 1 1 oc (n T ). Using the 
div-curl lemma, Lemma |3.10[ and Mu rat's lemma, Lemma |3.11| we show that 



3.15 



Thus we have established that 



RaxSax RS in ^ioc(nT), cf. Lemma 

(R-S) t -(c(u)(R + S)) x = 

holds weakly, cf. (3.51 1. By direct analysis of the scheme we infer that 

c{u) x — 2c(u)(R — S) weakly, 



cf. (3.49 1. Using the weak identity (c(u)u t ) x = (c(u)u x ) t we infer that u t = ^(R+S) 



holds weakly. To comp lete the argument, we derive an evolution equation for 



R 2 + S 2 , see Lemma 3.13 (cf. |T2l Lemma 11]) to conclude that u u — c{u)(c(u)u x 



weakly, and indeed that u is a weak solution. Our main result can be stated as 



follows (cf. Theorem 3.191: If uq and vq are such that i?( • ,0) and S( ■ , 0) are non- 
positive, and in i 3 (M) n L 1 (K), then the semi-discrete difference scheme produces 
a sequence that converges to a solution of ( |1.1| in the sense of Definition |1.1| 
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In Appendix |A] we show a higher integrability result, sec Lemma A.l which 
reads, here stated in the continuous case (cf. [12j Lemma 5]), as follows: 
H(R(-,0),S(-,Q)) G L 1 (R)ni 2 (M),thentt a; G L^ oc (Rx[0,T], c' (u)dx) for p G [2,3). 
The other results in this paper are independent of this, and the significance of the 
appendix is that it is suspected that such a regularity property could play a role in 
a uniqueness result. 



2. The difference scheme 



Our first aim is to construct an approximate solution of (1.51 based on a finite 
difference approximation of the spatial derivative. Rather than work on the full 
system of three equations, we derive approximate relations for the functions c(u) 
and c(it) in terms of R and S, thereby reducing the system to two equations. The 
temporal variable will not be discretized, and thus we will consider systems of 
ordinary differential equations indexed by the spatial lattice and depending on the 
lattice spacing. Subsequently we will show that as the lattice spacing decreases to 



zero, the system converges to the solution of (1.5 1. 

To avoid complicating the convergence analysis, we have chosen to restrict our 
attention to a semi-discrete difference scheme. To turn the difference scheme into a 
fully discrete one, we can rely on a variety of different time-discretization techniques, 
see Section [4] for one example. 



We shall use (1.51 as a starting point for a difference scheme. For j € Z, define 
Xj = jAx and Xj+1/2 = Xj + \Ax where Ax > is the lattice spacing. Let Ij 
denote the interval [2^-1/2, 2^+1/2)- 

Given any function K : K — > M, we let the value of K at the point Xj be denoted 
by Kj, that is, Kj — K(xj). 

On the other hand, given any sequence {Kj}j e %, we can consider it as the 
sampling at lattice points Ax Z of the function K defined by 

K{x) =Y d K j l Ij {x). (2.1) 

Here 1/ denotes the characteristic function of the set I. Clearly, if values {Kj} 
are computed from some difference scheme, we consider the function ( |2.1[ ) as the 
approximation of the true solution. 
It is easy to prove the inequalities 

\\ k \\l^{K) - IWIl 2 (r) ' II^Hl 2 (k) - H-^lli 1 ^) • 

Introduce forward and backward differencing by 

D±K j = ±± x {K m -K j ) 

for any sequence {Kj} of real numbers. Let (R,S) — {(Rj, Sj)}j e % satisfy the 
(infinite) system of ordinary differential equations 

R'jit) - c j+l/2 {t)D + Rj{t) = Cj(t) (R 2 3 (t) - S](t)) , (2.2) 
S'j(t) + c,.. 1/J (t)D.S j (() = -~ Cj (t) (R 2 j(t) - S](t)) , (2.3) 

for jeZ. The functions Cj ±1/2 and Cj are specified as follows. First set 

F(u) = [ 2c(v) dv. (2.4) 
J 
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Since c(u) > 0, we have F'(u) > 0, and F is therefore one-to-one. We start by 
defining Fj_ 1/2 by 

lim Fj_ 1/2 = 0, ^ 1 
^— or F j+1/2 = AxJ2(^-Si). (2-5) 

L>+*j-l/2 - Rj - 6j J ,=.00 

Then we can define Uj +1 / 2 by 

u j+ i /2 = (F- 1 ) (F j+1/2 ) , jeZ. (2.6) 

Note that this implies 

Rj - Sj = D+f 1 ^-.!^) = 2c (ut) D+Uj-i/2, 
for some value u" 1 " between Uj_i/ 2 and Uj+i/a- Therefore 

D+ »^ = fwy- (2J) 

Set 

Cj-l/2 = c(%_ 1/2 ), (2-8) 

and note that for some Uj between Uj_i/ 2 and Uj + i/ 2 we have 

c' (Vh 

#+c,-i/ 2 - c' («+) D +Wj -_ 1/2 = — V^p- (i?j - S,-) . (2.9) 

2C (Wj J 



So if we define 



c' («+) 



(2.10) 



J 4c («+)' 
we have that 

D + c j _ 1/2 = 2c j {R j -S j ). (2.11) 

Thus we have defined the functions Cj±i/ 2 — Cj±±/ 2 (R, S) and Cj — Cj(R, S). 
We will work with uq € if 1 (IR) and vq e L 2 (M). In this case we define 

^ju' (x)dx, r,,, ± 



Uo,j = uo(jAx), u' aj = —l u' Q {x)dx, v j = — j v (x)dx. (2.12) 



The initial values for (2.2 1 and (2.3 1 are 

Rj(G) =v 0J +c(u jfi )u' j, and Sj(Q) = v j - c(u j>0 ) u' j , (2.13) 
for j G Z. Extend the initial data {(Rj(0), Sj{Q))} je z by cf. ( pEl) , 

j 

S 0A x(x) = Sax(x,0) = J2S 3 (0)1i 3 (x). (2.14) 

At this point it is convenient to state the following general lemma. 
Lemma 2.1. Lei ip be a function in L 2 (M), and set 

( Pj=£^[ f(x)dx, (fiAx(x) = ^(fjlijix). 

X Jl i j 

Then 

\\tf - VAx\\ L 2 {M) -> 

as Ax — > 0. 
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Proof. For general functions <p,ip in L 2 (M) we have 

J (tpAx(x) - <PAx(x)) 2 dx = ^2 J J (^( Z ) ~ V{ z )) dz^ dx 

^E 



//, Ax J h 

- •' I 3 

(ip(z) - ip(z)) 2 dz. 



(i/j(z) - ip{z)Y dzdx (2.15) 



Thus 



\W~ <P&x\\ 2 < \\i> - <f\\ 2 + II^Ax - <P&x\\ 2 + IIV' - TpAx\\ 2 
< 2 \\ljj - lf\\ 2 + \\lp - %jjAx\\ 2 



(2.16) 



which shows that we, without loss of generality, can assume that tp is a smooth 
function with compact support, say supp(y) C [—M,M\. We find, similarly to the 
derivation of ( 2.15[ ), that 



((p(x) - (fAx(x)y 



dx = E^ (~Ex Jj ~~ tp{z))dz} dx 



^E 

3 



< 2 



1 

t. Ax 

M+l ( 



(ip(x) — <p(z)) dzdx 



(2.17) 



-M-i \ 2Ax J_ 



(tp(x) - ip(x - y)) dy\dx. 



Since tp is uniformly continuous, we can find <5 > such that \y\ < 6 implies 



\ip(x) -ip(x-y)\ < 



4(M + 1)' 



x e 



By choosing Ax < S we find that 



(ip(x) - (fiAx(x)) dx < e, 



which concludes the proof. 
This lemma implies 

|| Ro — -Ro,Az|| L 2 



^ as Ax -> 0, 
So - S ,Ax\\ L 2tm -> as Ax ^ 0. 



(2.18) 

(2.19) 

□ 

(2.20) 



Hypothesis 2.2. Consider u E W h3+q (R) (1 W 1 ' 1 (R) and v € L 3+9 (M) n L X (R) 
for some q > 0, and let Ro and Sq be defined by (1.6 1. Then we assume that Ro < 
and Sq < almost everywhere. 

This assumption implies that also Rj(0) and Sj(0) are nonpositive for all j. 
Furthermore, by interpolation, we have that uo G H /1,P (IR) and vo £ L P (R) for any 

pe [1,3 + ?]. 



Lemma 2.3. Assume Hypothesis 2.2 Then the system ( 2 .2 1 — ( 2 . 3 1 of ordinary 
differential equations with initial data (2.131 has a unique C l solution {i?j(f)} gZ 
and {Sj(t)}j eZ for all t > 0. 



8 



HOLDEN, KARLSEN, AND RISEBRO 



Proof. We use the notation R(t) — {Rj(t)} . gZ and S(t) — {Sj{t)} and write 
flZ2| and (Q, as 

i2j(t)=*?(JZ,S), 
S}(i)=*f(ie,S). 

Viewing this as an ordinary differential equation in ^ 1 (Z) x ^(Z), where £ 1 (Z) 
denotes the set of absolutely summable sequences with norm 

IMIli(R) = Ax^lVjl, 

i 

it will have a unique differentiate solution (R{t),S{t)) if 5) = , *f )} j& 
is locally Lipschitz continuous. This solution will be defined for t in some interval 
[0,t*) where t* is a "blow-up" time, i.e., 



If one can show that ||i?(i)|| i i^ R ^-|-||S'(t)|| i ij R -j < oo for all t > 0, then a continuously 
diffcrcntiable solution exists for all time. 
Now we claim that 



#(R, S) - tf(i?,S) 



< c 



R-R 



LH 



S-S 



LH 



where C is a constant that depends on || <S) ||xi(R) 5 1 1 S') Hi, 1 (m) 
shall show this for ^5> R , the arguments for fy s are identical. 
To show Lipschitz continuity we start by recalling (cf. (2.5)) 

j 

F j+ i/2 = F j+1/2 {R, S) = Ax & -Si). 



(2-21) 
and Ax. We 



Then 

Fj+i/2 (R, S) - Fj + i/ 2 (R, S) 



< 



Ax (\Ri-R. 



S{ S i 



< 



R-R 



LH 



s-s 



LH 



and therefore (writing Fj = Fj +1 / 2 (R, S)) 



F-F 




< 


R-R 


+ 


S-S 




L°°( 


R) 




LHR) 





LH 



Next we find (cf. ( pO) |) using ( |2~6] ) , pT8| ), and ( |2.11[ ) that 
~c J {R]-S 2 J )= l -D + c 3 _ 1/2 {R, J +S ] ) 



(cj+i/2 - c 3 _ 1/2 )(Rj + Sj) 



= ^(c(F ]+1/2 ) - C (F^ 1/2 ))(i?, + SO, 



abbreviating c(F J±1/2 ) = c((F 1 )(F j±1/2 )). Thus 



c^-S^-^R] 



^' 2 



(2.22) 



(2.23) 



i((c(J' i+1/2 )-c(F i+1/2 ))(iJ: i + 5 i ) 

- (0(^-1/2) -c^-.^^ + S,-) 
+ (c(F i+1/2 ) - cCFj.i/a))^ - ^ + - 4V 
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Now 



S) - 3f (£, S) < c (F J+1/2 ) - c(F j+1/2 ) \D + R 3 



1 



/)+ [Rj - Rj^ 



c(F j+1/2 ) - c(F J+1/2 ) (l^l + l^l) 

1 c(Fj_y 2 ) - c(Fj_ 1/2 ) (\Rj\ + \Sj\) 



c(F j+1/2 ) - c(Fj_ 1/2 ) ( R j - R 3 



S j S j 



< 



c 

Ax 



Fj+x/2 - F j+1/2 



C_ 

Ax 
C 



Ri — R'i 



\Rj\ + \R j+i \ 
Rj+i — Rj+i 



2Ax 

C 

2Ax 



F 



i+1/2 



F. 



i+1/2 



\Ri 



F 



j-l/2 - F j-l/2 



\Ri 



15,1) 



C 



Rj - Rj 



since c is Lipschitz continuous functions of Fj + i/ 2 . Multiplying the above by Ax 
and summing over j, we see that 



V R (R, S) - ^ R (R, S) 



< C 



where we have used (2.221 to find 



C = C AxA\R\\ 



\S\\ 



R-R 



R 



S~S 



s 



Therefore (2.211 holds, and we have established that {Rj(t)}j eZ and {<%(£)} - gZ 
exist for t < t * (fo r any initial data). If the initial data are nonpositive and in 



Lemma 



3.4 concludes the proof. 



□ 



Remark 2.4. The existence of global solutions of the system (2.2 1— ( 2 . 3 > with ini- 
tial data (2.12), that is, the fact that t* = oo, follows only after the estimate in 
Lemma 3.4 i.e., the inequality (3.131. Thus the results up to Lemma 3.4 are first 
valid for all times less than t* , and only after Lemma |3.4| we can infer that t* = oo. 
To simplify th e no tation, we state all these result for all t. 

For Lemma 2.3 we only require that uq € W 1,:L (R) and vo S i 1 (M). 

3. Convergence analysis 
Now let / be a sufficiently smooth function, and observe that 

/ (R j+1 ) = f (Rj) + f (Ri) (Rj+i - Rj) + \f"(rj) (R j+1 - Rj) 2 , 
where r 3 is between Rj+i and Rj. This can be rewritten 

D+f (Rj) = /' (Rj) D+R 3 + ^/" ( rj ) (D+Rj) 2 . (3.1) 
Furthermore, we have for any quantity /j, 

F>+ (cj-i/a/i) = Cj +x/2 D + fj + fjD+Cj_ 1/2 

= c j+1/2 D + fj + 2fjCj (Rj - Sj) . (3.2) 
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Similarly to (3.1 1 we have for a sufficiently smooth function g 

D_g (Sj) = g' (Sj) D_S 3 - ^-g" ( Sj ) (D_S 3 ) 2 , 

where Sj is between Sj and Sj—i. We also have 

D- (c j+1/2 gj) = Cj-yiD-gj + gjD+Cj-^ 



7-1/2-D-ffj + tgjCj (Rj - s^ 



(3.3) 



(3.4) 



Next define (cf. ( pTTTj ) ) 



ff(i?, S) = - (/' (R) - f (S)) (R 2 S 2 ) - (/ (R) - f (S)) (R-S). 
We shall use the following lemma repeatedly. 

Lemma 3.1. Let f,g € C 2 (R). Consider sequences {Rj}j^j,, {Sj}j£% satisfying 
Kty-Kty. For t > 0, there holds 



d f(Rj)-D + ( Cj _ 1/2 f(Rj)) + ^c j+1/2 f" ( rj ) (D + RjY 



dt 



2cj ( \f\Rj) (R 2 - S 2 ) - f(Rj) (Rj Sj) 



j t 9{S 3 )+D- (c j+1/2 g(Sj)) + ^Cj_ 1/2 g" ( Sj ) (D.Sj) 2 
1 



= ^9'(Sj) (R 2 - S 2 ) - g(Sj) (Rj - Sj) 
In particular, 

l(f( Rj )+f(Sj))-D + ( Cj _ 1/2 f(Rj))+D_ (c j+1/2 f(Sj)) 

+ ^ (c j+1/2 f ( rj ) (D + Rj) 2 + Cj . 1/2 f" ( Sj ) (D^S,) 2 ) 
= 2djH (Rj, Sj) , 
where rj is between Rj and Rj+i, and Sj between Sj and Sj— i- 



(3.5) 



(3.6) 



(3.7) 



Proof. Multiplying (j^) by /j = f'(Rj), using ([31]) and (j^), we find that 
J t fi ~ D + ( Cj _y 2 fj) + ^c j+1/2 f" (rj) (D +Rj ) 2 



= Tc 3 Q/j (R 2 - S 2 ) fj (Rj Sj)^ , (3.8) 



where fj = f(Rj). Similarly, multiplying ( |2.3[ ) with g 1 - = g'(Sj) for some function 
g E C 2 (R), using ([33]) and ^5A§, we find that 

j t 9j + D- {cj+i/29j) + ^Cj_ 1/2 g" ( Sj ) (D^S 3 ) 2 



[R] - S 2 ) ~ 9i (Rj - Sj^j , (3.9) 



where gj — g(Sj). Choosing f — g and adding (3.8 1 and ( |3.9[ ) we conclude that 
d3~7] holds. □ 

This lemma has several useful consequences, the first of which is the following 
result. 
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Corollary 3.2. We have that 

rT 



AxJ2 {R 2 j + S 2 ;) (*) + f AxJ2^ (c j+1/2 (D + R 3 f + Cj _ 1/2 (D_S 3 ) 2 ) dt 
j j 

<AxJ2(R 2 j+S 2 )(0). (3.10) 

3 

In particular, we have 

a x ( R2 + s 2 ) (*) < Ax J2 ( R i + s j) (°)- ( 3 - u ) 



Proof. Appy Lemma 3.1 with f(K) = K . In this case we observe that H(R, S) 



0, and /" = 2. Therefore, Lemma [3 . 1 1 yields 

J t (R 2 + S 2 ) - D + (c J - 1/2 i? J 2 ) + D- (c j+1/2 S*) 

+ Ax (c j+1/2 (D + Rj) 2 + Cj _ 1/2 {D_Sjf) < 0. 

Multiplying with Ax, summing over j, and integrating in t finishes the proof of the 
corollary. □ 

The variational wave equation enjoys certain invariance properties in the (i?, S) 
variables. Indeed, if both are nonpositivc initially, they will remain so for all time. 
Furthermore, if in addition the initial data are bounded below by a (negative) 
constant, then the same constant bounds the solution for all time. See, e.g., |T51 
Thm. 3.1.6]. The approximate solution has the same properties, which is the result 
of the following lemma. 

Lemma 3.3. The following statements hold: 

(i) IfRj(0) < and Sj(0) < for all j, then Rj(t) < and Sj(t) < for all t > 
and all j . 

(ii) If —M < Rj(0) < and —M < Sj(0) < for some positive number M and for 
all j, then -M < R 3 (t) < and -M < Sj(t) < for all j and t > 0. 



Proof. To prove the first statement (z) choose f(K) = (max{0, K}) in (3.7 1; with 
this choice 

{0 if RS > 0, 

RS(R -S) if S < and R > 0, 
RS(S -R) if S > and R < 0. 
Hence H (R, S) < 0, furthermore f"(K) > 0, and thus using ( |3.7[ ) we find that 

((max-jO,^)}) 2 + (max{0,5j(t)}) 2 ) < 0, 

3 

since Rj(0) < and Sj(0) < for all j. Thus the first statement (i) of the lemma 
holds. 

To prove the second statement (ii) choose f(K) = (min {K + M, 0}) . Then we 
find that 

'0 if R > -M and S > -M, 

-2M(R - S) 2 if R < -M and S < —M, 

(R + M) (R — S)(S — M) if R <-M <S, 
(S + M) (R — S) (M — R) if S < -M < R, 
which implies that 

H(R,S) \ {(RiS)lR<MiS<M} < 0- (3-12) 



H(R,S) 
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Furthermore f"(K) > 0. Thus, if > Rj(0) > -M and > Sj(0) > -M, we 
observe from the first statement (z) that Rj and Sj remain negative. This implies, 
using (3.121, that H(Rj, Sj)(t) < 0. Hence it follows as before, using equation 
pj| , that 

( (min {0, Rj(t) + M}) 2 + (min {0, S 3 (t) + M}) 2 ) < 0. 



Thus the second statement (ii) of the lemma follows. □ 

In case Hypothesis |2.2| holds, we have the integrability estimate. 
Lemma 3.4. If Rj(0) < and Sj(0) < for all j, then 

Axj2(\m)\ p +\m\n<^J2(\ R ^\ p+ \ s M p ), (^) 



for any p > 1. In addition, if Hypothesis 2.2 holds, then for any p € (2, 3 + q) 

[ AxVc'(ut) \D+u^ 1/2 \ p+1 dt<C p , T , (3.14) 
Jo j 

where C Pi t is a constant depending on p and T (but not on Ax). 

Remark 3.5. This lemma finishes the proof of Lemma |2.3| namely the fact that 



t* = oo, cf. Remark 2.4 



Proof of Lemma \3A\ Choose f(K) — \K\ V , and observe that 

f(0) = 0, f(K) = f(-K) and f"{K) > 0. 
Now it is easy to see that 

H(R,S) = H{S,R) and H(-S, -R) = -H(R, S). 

Furthermore 

H(R,R) = H(R,-R) = 0. 

We also find that 

VH(R, S) ■ (1, 1) = \ (f"(R) - /"(£)) {R 2 - S 2 ) > 0, 
since /" is an even non-negative function. From this it follows that 



H(R, S) 



<0. 



\{(R,S)\R+S<0}- 

Hence, since So < and Rq < 0, by Lemma 3.3 also Rj(t) and Sj(t) are nonpositive 
for t > 0, and thus 

AxJ2(f(Rj(t)) + f(Sj(t))) <Ax^ (f(Ri(P)) + f(Sj(0))). 



For the proof of (3.14), we fix p £ (2, 3 + q), remember that Rj < and Sj < 0, 
and calculate 



l r 



H(R j ,S j )=^ p ( sign (R^lRjf- 1 - sign (S,-) {tij - 



ip-i 



l r 

2 

1 

2 



2{\R^~\S^) {Ri-S-) 

-p^r'-is.r 1 ) (i^-i 2 -i^i 

f2(|i? J i-i5,i)(ii? J r~i^n 
- P (\R J \^\s J \) 2 (\R J r i + \s J r i 
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21^11^1(1^1-1^1) l\R 



IP-2 



IP-2 



+ 2(\R,\~\S J \)(\R ] \ P -\S ] 












+ b{R 3 ,Sj) 


+ 2(\R j \-\S j \)(-\R j \\S j \ 


0* 




-w a ) + w 


' r p(\R,\-\S 3 \) 2 {\R/-U 


i^r 1 ) 




+ 2(\R J \-\S J \) 2 (\R J \ p - 1 + 






-( J) -2)(\R j \-\S j \) 2 {\R j 




i +\s j r i )+b(R j ,s j ) 



_ 1 

~ 2 

It is easy to see that &(i?, S) < for p > 2, and we also have the inequality 
for some positive constant K p depending on p. Hence 



H (R h Si) < Kp( P 2) \Rj - S.r 1 = -K p (p - 2)c («+) \D 



+ Uj-l/2\ 



I P+1 



By Hypothesis |2.2| we find that 
K p {p-2) 



uj) \D+Ui- 



+Uj-l/2\ 



IP+1 



< 



from which (3.141 follows. 



Extend the functions (Rj,Sj) to the full line, cf. ( |2.l| and ( |2.14| , by 
RAx(x,t) = ^R j {t)l Ij {x), and S A x(M) = ^ S,-(t) 1/^ (a:). 



Define Fa x by 



and then uax by 



FAx(x,t)= / (R Ax (x,t) - S A x(x,t)) dx, 



L 



2c(u) = FAx{x,t). 



AxY^(\RM\ p +\sm p )<c, 

3 

□ 

(3.15) 

(3.16) 
(3.17) 



Note that 



or 



D+F Ax (xj-.i/a, t) = flj - = Z?+F (^- 1/2 (t)) , 



2c(w) cfo = 



i-1/2 



«Ai(lj + l/2,*) 



'»A.(lj-l/2:*) 



2c(i>) dv. 



Now we have that lim-r^^oo Rax(x, t) = lim. E ^_ 00 Sax(x, t) — 0, and therefore 
lim.x—,-00 ua x {x 1 t) — 0. Hence we must have UAx{xj^i/2,i) = Uj_i/ 2 (t) for all j. 
It is convenient also to define the piecewise constant function 



u A x(x,t) = ^2u J _ 1/2 (t)l Ij _ 1/2 (x). 



(3.18) 



Now we can show the (local) uniform convergence of UAa 
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Lemma 3.6. Assume Hypothesis |2.2| Then there exists a function u G C(Ht) 
such that for any finite interval [a, b] we have 

lim UAi(i,t) = u(x,t) uniformly for (x,t) £ [a, 6] x [0, T]. 

A.r 



(3.19) 



\\RaA • ,*)|| l1(r) + ||£a*( • , t)|| L i (M) < ||i?Ax( • ,0)|| L i (R) + \\s Ax ( ■ , o)|| il(R) < c 

(3.20) 

where C is the constant in (3.191. This implies that F Ax is uniformly bounded, 
since 



Proof. From Hypothesis |2.2| we infer that 

||-Rax( • ,0)|| iP(R) + \\S Ax ( ■ ,0)\\ L p(r) < c, 

for both p = 1 and p = 2 for some constant C that is independent of Ax. From 
Lemma \3. 41 it follows that 



(RAx(y, t) - S Ax (y, *)) rfy 



< C. 



1^(2!, «) | = 

Next, we observe that 

F Ax (x,t) - F Ax (0,t) = f (R Ax (y,t) - S Ax (y,t)) dy. 
Jo 

Therefore, using Cauchy-Schwarz's inequality and ( |3.20| , we find 

\\F Ax (-,t)-F Ax (0,t)\\l 2(ab) = J (J (R Ax (y,t)~S Ax (y,t))dy^j dx 

< f 2\x\ f (R 2 Ax (y,t) + S 2 Ax (y 7 t))dydx 

J a JR 

< C 2 (a 2 + b 2 ) . 



Similarly, by using (3.161 we find that 



\\d x F Ax (-,t)\\ 2 L2{ab) < (ll Ax (x,t) - S Ax (x,t)\ dx 



< 2 



\R Ax (x,t)\ + \S Ax (x,t)\ )dx 



< 2(a + b)C 2 



using (3.11 1. Thus there is a constant C±, independent of t and Ax (but depending 
on a, b), such that 

Morrey's inequality now implies that for x and y in [a, b] we have that 

\F Ax (x,t) - F Ax (y,t)\ <C 2 \x-y\ 1/2 , (3.21) 

for some constant C 2 which is independent of Ax and t (but depending on a, b). 

Now note that using f(R) = R in ((OJ and g(S) = S in we find, cf. $t~ 
that 



R'j - D + { c j-i/2Rj) - -cj (Rj - SjY 
Sj + D- (c j+1/2 Sj) = -Cj (Rj - SjY 



Since D + Fj — Rj — Sj, 
d 



3-1 



dt 



Fj = Ax ]T (R'i - Si) = Cj _ 1/2 (Rj + Sj-!) . 
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Therefore we have that for < s < t < T 

d 



\\F Ax {-,t)-F Ax (^s)\\ L2W < f 

J S 



<c. 



F Ax {-,r) 



\Sax( ■ 



dr 



\lh 



\Rax(-,t) 



\LH 



dr 



< C 3 C \t - s\ , 



using (3.11 ). 

Since H i (a, b) CC C{a,b) C L 2 (a,b), we can use [HI Lemma 8] to deduce that 
for x and y in (o, b), we have that for any rj > 0, there is a finite > such that 

\F A x(x,t) - F Ax (x,s)\ < r)\\F Ax (- ,t) - F Ax ( -,s)\\ H1{ab) 

+ C v \\F Ax ( ■ ,t) - F Ax ( ■ ,s)\\ L 2( atb ) 

< r)2d + C V CC 3 \t-s\. 

For any e > we choose (x, t) and (y, s) in [a, b] x [0, T) and rj > such that 

C 2 \x-y\ 1/2 < J, ?72Ci < \ and then C n CC 3 \t - s\ < J. 

With this choice 

\F Ax {x,t) - F Ax {y,s)\ < \F Ax (x,t) - F Ax (y,t)\ + \F Ax {y,t) - F Ax (y,s)\ 
< e. 

Hence, the sequence {F Ax } Ax>0 is equicontinuous and uniformly bounded in [a, b] x 
[0, T], and by the Arzela-Ascoli theorem there exists a convergent subsequence 
(which we do not relabel). 

By the definition (3.171 of u Ax and the assumption on c, cf. ( |1.9[ ), we find that 

\F Axj (x,t) - F Axk (x,t)\ > d \u Axj (x,t) - u Axk (x,t) \ , 

for some constant C4 depending only on the function c. This shows that {u Axj } is 
Cauchy and thus uniformly convergent on compacts [a, b] x [0, T]. □ 

Remark 3.7. For this lemma to hold it is sufficent to assume that Rq and Sq (and 
therefore i?Aa:(0), Sa^O)) are nonpositive, and in L 1 n L 2 . 

Note that F Ax is linear in the interval Ij as R Ax and S Ax are constant there. 
By definition we have that 

du Ax 1 



dF Ax 2c(u Ax ) 



> 0. 



This means that u Ax (x, t) is monotone in the interval Ij, and we have that u Ax (Xj±i 
Uj±i/2(t). To simplify the subsequent calculations we introduce 



/2: 



and 



= 0j u j-i/2 + (1 - 0j)u j+1 /2, 9j €[0,1], j el, 

u Ax (x,t) = y^uj(t)i I:i (x). 



Then for any fixed x and t, 



lim u Ax {x,t) — u{x,t). 



(3.22) 



This is so since if x € Ij, there is a yj € Ij such that u Ax (x, t) — u Ax {yj, t) by the 
monotonicity of u Ax . Now let {Aa;^} and {Aa;^} be two sequences tending to zero. 
Fixing x, we can find sequences {yk} and {yi} such that y^ — > x and yi — * x and 

u Axk {x,t) = u Axk (y k ,t) and u Axi (x, t) = u Axr (y e , , t) . 
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Hence 

\uAx k (x,t) - UAx,(x,t)\ < \u Axk (yk,t) - u Axe (yk,t)\ + \uAx t (Vk,t) - u Axe (ye,t)\ . 

Both terms on the right vanish as k and I become large since u& x is uniformly 
continuous. Hence for any choice of {6j}, ( |3.22 1 holds. In particular, this implies 
the pointwise convergence 

lim ^""^ Uj± 1 / 2 (t)li j (x) — u(x,t) and 

uniformly on compacts. 

Next, we collect (in three lemmas) some well-known results related to weak 
convergence. Throughout the paper we use overbars to denote weak limits. 

Lemma 3.8 (|2J. Let O be a bounded open subset ofM. M , with M > I. 
Let {v n } n>1 be a sequence of measurable functions on O for which 

sup / $(\v n (y)\) dy < oo, 

n>l JO 

for some given continuous function $ : [0, oo) — > [0, oo). Then along a subsequence 
as n — > oo 

for all continuous functions g : K — » K satisfying 

iim S = °- 

|v| — *• oo $(H) 

Lei g: M — » (—00,00] be a lower semicontinuous convex function and {t>n}„>i a 
sequence of measurable functions on O, for which 

v n v in L l (0), g(v n ) £ i 1 (0) for each n, g(v n ) g(v) in L 1 (0). 

Then 



g(v) < g(v) a.e. on O. 

Moreover, g(v) € i 1 (0) and 

/ g(v) dy < liminf / g{v„) dy. 
Jo n ^°° Jo 

Lf in addition, g is strictly convex on an open interval (a, b) C M and 

g(v) = g{v) a.e. on O, 

then, passing to a subsequence if necessary, 

v n {y) -> v(y) for a.e. y e {y E O \ v(y) £ (0,6)}. 

Let X be a Banach space and denote by X* its dual. The space X* equipped with 
the weak-* topology is denoted by X^ oak , while X equipped with the weak topology 
is denoted by X wca k- According to the Banach- Alaoglu theorem, any bounded ball 
in X* is <t(X*, X )-compact. If X separable, then the weak-* topology is metrizable 
on bounded sets in X*, and thus one can consider the metric space C ([0, T];X* cak ) 
of functions v : [0, T) — > X* that are continuous with respect to the weak topology. 
We have v n -> v in C ([0, T];X* cak ) if (v n (t),4>)x*,x -> («(*), <^>x*,x uniformly 
with respect to t, for any e A. The following theorem is a straightforward 
consequence of the abstract Arzela-Ascoli theorem: 
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Lemma 3.9 ([2 ). Let X be a separable Banach space, and suppose v n : [0,T] — * 
X* , n = 1, 2, . . . , is a sequence of measurable functions such that 

II v "IIl°=([o,t] ; x*) - 
for some constant C independent of n. Suppose the sequence 
[0,T] 3 t« (v n (t),$) x *,x, n = l,2,..., 

is equi- continuous for every $ that belongs to a dense subset of X . Then v n belongs 
to C ([0, T];X* eak ) for every n = 1, 2, . . . , and there exists av£ C ([0, T];X*^) 
such that along a subsequence asmco 

v n ->v in C([Q,T];X* eak ). 

Lemma 3.10 (Div-curl lemma 7 ). Let QcR 2 be a bounded domain. Suppose 

vl^v 1 , v 2 ^v 2 , 

w\ —^w 1 , w 2 —^w 2 , 

in L 2 (Q) as e — > 0. Suppose also that the two sequences {div (yl, i» 2 )} E >o and 
{curl (w',Wg)} £> o lie in a common compact subset of H^(Q) , where div (yl,v 2 ) = 
d Xl v\ + d X2 v 2 and curl (wl,w 2 ) = d Xl w 2 — d X2 w\. Then along a subsequence 

{yl,vl) ■ (wl,w 2 ) -> (V,v 2 ) ■ (wJ 1 ,^ 2 ) m V'(Q) ase^O. 

Lemma 3.11 (Murat's lemma |7]). Suppose that {£ £ } e>Q is bounded in W^'°° (Ht) ■ 
Suppose also that C e = C\+C 2 , where {£e} £>0 lies in a compact subset of H^(Ht) 
and {£ 2 } £>0 lies in a bounded subset of A4\ oc (Ht) ■ Then {C e } £>0 lies in a compact 
subset o/ff^pT). 



Accor di ng to Hypothesis ^ R ,S £ L 1 (R) n L P (K) (with p > 3). In view 
of fl3.10| , ( |3.13| ) and Lemma |3.8| there exist R, S £ L°°(0, T; L q (R)), q £ [l,p], 



R 2 , S 2 £ L°°(0, T; L r (R)), r £ [l,p/2], such that along a subsequence as Ax -> 

R Ax ^ Rin L oc (0,T;L 2 (R)), Sax S in L°°(0, T; L 2 (]R)), 

(3.24) 

#L ^ ^ 2 in i°°(0,T;£ r (]R)), 5^ A S* 2 in L°°(0, T; L r (R)). 

As a matter of fact, we can assume that for any function / £ C 1 (M), with 

\f(z)\ < C (l + |z| 2 ) and |/'(z)| < C (1 + |«|), (3.25) 

the following statements hold 

/(#A*) A TP), /(Sa*) A 7(S) inL°°(0,T;LP/ 2 (K)), (3.26) 

where the same subsequence of Ax — + applies to any / from the specified class. 
Clearly, we can also assume that as Aa; — > 

5(^Ax) (Rax - S Ax f -± c(u)(R - S) 2 = c(«)p - S) 2 inLP/ 2 (n T ), (3.27) 



by equation (3.23). From Lemmas 3.1 and 3.4 it is not difficult to deduce that the 
functions 

t» [ f(R A x)$dx, tn f f(S A x)<S>dx (3.28) 
are equi-continuous on [0, T] for every $ £ Cq°(R). In addition, f(RAx) an d f(SAx) 



are bounded in L°°(0, T; L' r (R)), independently of Ax. Now we apply Lemma 3.9 
with X* = L r (R), X = L r ' (K), and r' = r/(r- 1). Since Cfi°(R) is dense in 27 (M), 
we can thus assume that f(R),f(S) £ C ([0, T]; L^akW) and 



/(i?Ax) - /(i?), /(^Ax) - in C([0,T];L; cak (]R)). (3.29) 
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Of course, when f(z) = z, we can assume R,S £ C ([0, T]; i^ cak (E)) and 

Ra x ^R, S Ax ^SmC([0,T];L 2 weak (R)). (3.30) 



Lemma 3.12. Assume Hypothesis\2.2\ Then we have, cf. (1.8 1, 



R t -(c(u)R) x = -c(u) (R-S) 2 , 



and 



S t +(c(u)S) x = -c(u) (i?-S) 2 , 
in the sense of distributions onRx [0,T), i.e., for any (p £ C^°(K x [0,T)), 

(Ript — (c(u)R) (fx) dx dt + / Ro(x)(p(x, 0) dx 



(3.31) 
(3.32) 



c(it) (R — S) 2 ip dx dt 



and 



[Sipt + (c(u)S) ip x ) dx dt + / So(x)(p(x 7 0) dx 



[[ c(u)(R- S) 2 dxdt. 
J Ju T 



Proof. Fix ip € Cg°(K x [0, T)). The equation with f(R) = R reads 



dt 



(3.33) 



Set 



Ax 



ip(y,t) dy. 



Next multiply the equation (3.331 with ipj, sum over j, do a partial summation, 
integrate over i, to end up with 



{RaxVi - caxRax^x) dxdt + / Ro,Ax(x)ip(x, 0) cfe (3.34) 

R 

cax(Rax - S Ax ) 2 tpdxdt 

+ / y^Cj-Rj / (D-ifj — (fix) dxdt 
Jo j J, 

where we have defined the functions c Ax and cai by 

c Ax (x,t) = ^2c j _ 1/2 (t)l Ij _ 1/2 (x) and c Ax (x, t) = c j {t)l Ij (x). 



By (3.231 



Now 



cai — ► c(it), cai — > c(m) uniformly on supp (ip). 



(3.35) 



l^-ViC*) - ¥>*0M)l < ll^|lioo (Kx[ o, T ]) Ax, x e Ij. 
In view of this and ( |3.13 1, the last term in (3.341 is bounded as follows: 

/ y^Cj-Rj / (D-ifj — (fx) dxdt 
Jo , •// 



< CAx -> 0. 



Furthermore, in view of (2.201, as Ax — > 



^o,Aa;( :z; ) l i 9 ( :z 'i 0) ^ — > / Ro(x)f(x,0) dx 
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Hence, sending Ax — > in (3.34) yields (3.311. The evolution equation (3.32) 
for S is proved in the same way. □ 

We can also prove a generalization of the previous lemma. 

and let f G C 2 (R) be a convex function 



Lemma 3.13. Assume Hypothesis 
satisfying (3.25). Then 



2.2 



f(R) f - (c(u)f(R)) x < 2~c(u) Q/'(i?)(i? 2 - S 2 ) - f(R)(R - S)) , (3.36) 

f(S) t + (c(u)J(Sj) x < -2c(«) Q/'(S)(i? 2 -S 2 ) - f(S)(R-S)) , (3.37) 
in the sense of distributions onlx [0, T), i.e., for any tp G C§°(IR x [0, T)), <p > 0, 



c(«)/(ii) p a da;dt + / /(i? o (^))^,0)da: 



2c(u) { -f'(R)(R 2 - S 2 ) - f(R)(R - S) ) (pdxdt 



and 



f(S)<p t + (c(u)f{S))ip x ) dxdt + f(S (x))<p(x,0) dx 



2c(u) ( ^f'{S){R? - S 2 ) - f(S)(R - S) ) tpdxdt. 



Proof. Similar to the proof of Lemma 3.12 starting from (3.8 1 and (3.9 1. 

The weak limits R 2 , S 2 satisfy the initial data in a strong sense: 
Lemma 3.14. Assume Hypothesis\2.2\ Then 



a 



lim 



lim / R 2 dx = lim / R dx — / R?.dx 



S 2 dx — lim I S dx 



(3.38) 



S dx. 



Proof. Since R,S G C ([0, T]; L 2 wcak (R)), it follows from ( pOl) ), fl3~32| that 
i?(-,t) ->> i? , in L 2 (M) as £ 0+. 



From this, (3.241, and Lemma 3.8 we conclude that 



/ Rq dx < lim inf / R 2 dx, / S 2 dx < lim inf / S 2 dx. 



(3.39) 



On the other hand, ( |3.24[ ) says that R A x(-,t ) R{; t), Sa*(-,*) in L 2 (M) 

for a.e. i > 0, and thereby, using also (3.101 and ( 2.20[ ), 

(if + S' 2 ) (t, x) < 

Since R 2 , S 2 G C ([0, T]; L^, eak (M)) (with r > 1), one can prove that this inequality 
actually holds for alH > 0. Combining (|3.3Qb and (|3.40| yields (|3.38|. □ 



(r 2 + S 2 ) (i, cc) dx < { (R 2 + S$) dx. (3.40) 



Lemma 3.15. Assume Hypothesis 2.2 
satisfying \f(z)\ , \g(z)\ < C\z\. Then 



and let and let f, g G C 2 (R) be functions 



f(RAx)g(SAx) ~ * f(R)9(S) i n the distributional sense on M. x (0,T). (3.41) 
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Proof. We will show that the sequences 

{dtf(RAx) - d x (c(uax) f(RAx))} Ax>0 , 

{d t g(S A x) + dx (c (u Ax ) g(S Ax ))} Ax>0 

are compact in H{£(R x (0,T)). 

Introducing the distribution C Ax — dtf(RAx) — d x (c(u Ax ) /(Rax)), we find 

(£ax,<p) 



- r A i [ \ f'{RAx) {Rlx ~ S 2 Ax ) - J(R Ax )(Rax - S Aa 



c f 



/ y2 c jf( R j) (D-tpj - ip x ) dxdt, 
Jo „■ Jn 



ip dxdt 
(3.42) 

(3.43) 



for <j) e C °°(Mx(0, T)), w here C f (x, t) is a function that is bounded in L x (Mx(0, T)) 
independently of Ax, cf. (3.34). The last term above is bounded by 

r T ^ 

L/j — <-PxW ro m ~. dt. 





PAx ■= C\\Ro\\ L 2{R) / ||2J- D -^' 1j 

Jo 



2.1 



Since D-ipjlj^ is a piecewise constant approximation to ip x , by Lemma 
p Ax tends to zero as Ax — ► 0. Thus we infer that 

\(£a x ,<p)\ < C|MIl°°(Rx(0,T)) +PAx, 

where p Ax tends to zero with Ax. Thus (3.431 is in a compact subset of ff,~^(K x 
(0, T)), while (3.42 1 is in a bounded subset of the locally bounded Radon measures. 
Hence Murat's lemma implies that C Ax is compact in Hr~(M. x (0,T)). By analo- 
gous arguments, {d t g(S Ax ) + d x (c (u Ax ) g(S Ax ))} is compact in (M x (0, T)). 
Now by the div-curl lemma on the sequences 

{g(S Ax ),c (u Ax ) g(S Ax )} Ax>0 and {c(u Ax ) f(R Ax ), f(R Ax )} Ax>0 , 

we see that 



2c{u Ax ) f{R Ax )g{S Ax ) — > 2c(u)f(R)g(S) in the distributional sense, (3.44) 
which, due to ( 3.44[ > and Lemma 3.6 concludes the proof of (3.41). □ 

An immediate consequence of the previous lemma is the following result . 

Corollary 3.16. There holds 

(R-S) 2 = R 2 -2RS+S 2 a.e. in R x (0,T). (3.45) 

Proof. Since we can assume that R Ax S Ax RS in L 1 (M x (0,T)), it follows from 
Lemma 13.151 that 



x(0,T) 



RS ip dxdt — 



x(0,T) 



R Sip dxdt, Wip G C °°(R x (0,T)), 



from which we infer that RS = RS a.e.; Hence (3.45) follows. 
Lemma 3.17. Assume Hypothesis\2.2\ Then 



R 2 = R 2 and S 2 — S 2 , for a.e. (x,t) G Ht- 

Consequently, as Ax — > 0, 

-Raz - * -R) — > S m L 2 oc (Ht) and almost everywhere in Ht- 

Proof. Using Lemmas |3.12||3.14| and Corollary |3.45| we can argue exactly as in, 
e.g., Zhang and Zheng [2], to arrive at (3.461 and (3.471. □ 



□ 



(3.46) 



(3.47) 
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Lemma 3.18. Assume Hypothesis 2.2 Then u is a weak solution of (l.ll, i.e. 



d u 
__ c(w )_ (c(u) 



dx 



= 



weakly in LTy in the sense that 

(ut<p t - c(u) x (c(u)tp) x ) dxdt = 



(3.48) 



for all test functions (p G Cq°(Ht)- Here u t and c(u) x are given by (3.531 and 
(3.491, respectively. 



Proof. We claim that 
To this end let 



c(u) x = 2c(u)(R — S), weakly. 

CAx = ^Cj-_i/ 2 lzp 



(3.49) 
(3.50) 



and compute 



dx 



CAx, V 



cax<Px dxdt 
-J ^ J Cj-x/2<p x dxdt 

y2 C J-l/2 D + l P ( x J-l/2, t) Ax dt 

Jo j 

f T 

/ 2_. (£>+Cj_i/ 2 ) (fi (afj-_i/ 2 ,t) Axdt 
Jo . 

[ AxJ2^ *) [2c, (Rj - S )} dt. 

Jo 



By sending Ax to zero in this equality, and using (3.23), our claim (3.49) follows. 
From Lemma [3.12| we find that 

[R — S) — (c(u)(R+ S)^j x = 0, in the sense of distributions. (3.51) 

Observe that for a function u that is at least one time differentiable we have that 
(c(u)ut) x = (c(u)u x ) t holds in the distributional sense. Specifically, we have 



(c(u)ut)(p x dxdt = // c(u)tp x ut dxdt = // ((c(u)(f) x — c'{u)u x <p\ut dxdt 

Tt «* 

[(c(u)<p)t — c'(u)ut(p)u x dxdt 

(c(u)u x )(ft dxdt. (3.52) 
Thus we see that this can be rewritten 

~(c(u) (2u t -(R + S)))=0 
in the sense of distributions. Hence 

u t =*(R + S), 



(3.53) 



since 



lim u(x,t)= lim (R(x, t) + S(x, t)) = 0. 



x — * — oo 



X X 
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Set 

R £ (x,t) = f R(y,t)f(x-y)dy, 
Jm 

where j e is a standard mollifier. Then 

Tf t - (c(u)R £ ) x = -c(u)(R~S) 2 * f + r e , 
where by the DiPerna-Lions folklore lemma 

r s = (c(u)R) x * f - (c(u)R% 
tends to zero in L\ oc {YIt). This in turn implies that 

R e t - c(u)R e x = -c(u)(R-S) 2 * f + 2c(u)(R - S)R £ 



= -~c{u) ({R - S) 2 - 2R(R - S) 
Similarly, with S £ = S * j £ , we get 
c(u)S £ x 



c(u)[(R-S) 2 + 2S r(R-S) 



(3.54) 



(3.55) 



where s £ tends to zero in L1 oc (Ht). Adding (3.54) and (3.551 and sending e to zero, 
we get, after using Lemma |3.17| 



d_l 
dt2 



! U S i r(u) ^\( R ~ S ) = 5 ( u ) {{R- s) z - (R- S) 2 ) =0, 



which, using ( |3.53[ ), (3.491 and Lemma [3. 15 can be rewritten as ( |3.48 1. 
We collect some of our results in the following theorem. 



□ 



Theorem 3.19. Assume (1.9 1 and Hypothesis 2.2 Then the semi-discrete differ- 
ence scheme defined by (2.2 1 — ( 2.10 ) produces a sequence that converges to a weak 
solution of (1.1 1. 

The same conclusion holds if Hypothesis 2.2 is replaced by the assumption that 
uq G V4 7l,1 (K) and vq 6 i 1 (IR) and that Rq and So take values in [— M, 0] for some 
positive constant M . 

Proof. Observe that in the case when Rq and Sq take values in [— M, 0], Lemma 



3.3 shows that Rj(t),Sj(t) € [— M, 0] for all t. Thus R Ax ,S A x € L°° n L 1 , and by 



interpolation we see that Hypothesis |2.2| is satisfied 



□ 



Remark 3.20. The equation (2.7) is awkward to use in practice, since one must 
compute the inverse of F at each node. In order to circumvent this we may redefine 
Cj slightly. Let Rj and Sj be defined as before, but let Uj_ 1 / 2 be defined by 



D+Uj-l/2 



R i 



7 x- ( 3 - 56 ) 

c (U 3 -l/2) + C{Uj+l/2) 

This is also a nonlinear equation to solve for Uj+1/2, but solving this is likely to be 
easier than inverting F. However, if we accept a certain imbalance, we can define 



/2 by 



D + u J-l/2 = 



2c(u j _ 1/2 )' 



(3.57) 



In order to get our approach to work, it is essential that (2.111 holds. Therefore, 
; this is the case. Since c is contin 

D + c(uj-_i/ 2 ) = d (uj) D + Uj-i/2, 

If Uj + 

c'(%) 



we shall define Cj so that this is the case. Since c is continuous, we have that 



for some Uj between Uj_i/ 2 and Uj +1 / 2 - If Uj +1 / 2 is defined by (3.561 then we set 



2 (c(uj-_i/ 2 ) + c(uj+i/ 2 )) ' 



(3.58) 
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(3.59) 



while if Wj+i/2 is defined by (3.571 we set 

1 4C ( W J-l/2) ' 

In both cases (2.111 holds. Therefore, the schemes defined by (2.2 1, (2.2 1, and 
(3.561, (3.58) or (3.571, (3.591 all produce sequences converging to weak solutions 
ofTLl). 



4. Numerical examples 



The semi-discrete scheme defined by (2.2 1— ( 2.12 1 is rather involved, in particular 
the computation of the quantity Cj. For actual computations one needs to make 
a further discretization of the time variation. We have considered the following 
versions of the semi-discrete scheme defined by (2.2 )-( 2.11|. These schemes all use 

(4.1) 



an explicit discretization of (|2.2|)-(|2.3|) 

c 



1+1/2 D + R] = c] ((i??) 2 - (S7) 2 ) 



where 
and 



C j-l/2 - C (""-1/2) 



D + K j 



Furthermore, since we wish something like ( |2.11 1 to hold, 



c" 



2{R^ - SJ) ' (4 ' 2) 

The difference between these schemes consists in the way the "cofficients" c™ 
are defined. 

(1) Integration in time. We update u™_ 1 , 2 by considering a discrete version 
of ( |3~53| , 

J + l/2 - "i + 1/2 

+ y [R] + + + Rjtt + s " + s ?+i + s ? +1 + s ?ti) ■ ( 4 - 3 ) 

We use this scheme since u is discretized on a grid that is staggered with 
respect to that of R and S. 

(2) Integration in space. Knowing u™ N _ 1 ^ 2 for some large N, we can set 



U ]+l/2 — U j-l/2 + 



cm on 



C K-l/ 2 ) + C ( U "+l/2)' 

In this section we describe two examples. 

Consider first the case where the function c is given by 

2 

cOu) — — (-7T + arctan(u)) , 

7T 

and the initial data are given by 

R(x,0) = -2e- (x - 5 ) 2 , S(x,0) = 2e' (x+5 ^ . 

In this case 



(4.4) 

(4.5) 
(4.6) 



u(x, 0) 



x R(y,0) - S(y,0) 1 . 

— — dy, u t (x,0) = -(R(x,0) + S(x,0 ). 

-oo 2c( u (y ; 0)) 2 



In Figure [T] we show the computed solution u (top) and R and S (bottom). The 
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x 



Figure 1 . The scheme ( |4.4| with the initial data ( |4.6| and c given 
by (4.5). The u variable (top), the R variable (middle) and the S 
variable (bottom) as functions of x and t. 



discrete difference scheme can be studied numerically also in cases not covered by 
the convergence results in this paper. We have included an example of that type 
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here. Define the function c by 



c{u) = ^Jacos 2 (u) + /3sin 2 (u), a = 1.5, = 0.5. (4.7) 
When testing, we take the initial data from [3J, and use 

u(x,0) = J + e~ x2 , u t (a! > 0)=-c(u(a:,0))-^-u(a! > 0). (4.8) 
4 aa; 

In order for the two schemes to be compatible, we have defined u^/j by 

o 

D+u J- 1/2 = c(u° +1/2 ) + c( u °_ 1/2 y 

even for the scheme using by fl3~3) . In Figure |] we show u for the two methods 
with initial data (4.8 1 using Ax = 30/256, and At = Ax. We remark that using 
At = Ax/M where M is a large integer, produced very similar results. 

Appendix A. Higher integrability properties 

In this appendix we prove a so-called higher integrability result. Briefly stated, 
we have that if Rq and So are nonpositive and in L 1 n L 2 , then d x u{ ■ , t) is in Lf oc 
for all p e [2,3). This is obvious if Rq and So are in L 3 , and the significance of 
this section is that the d x u is more integrable than is to be expected. The reason 
for including this is that we suspect that such a property will play a role in a (yet 
unknown) uniqueness result. 

Throughout the appendix we assume that Rq and So are nonpositive and in 
£ 1 (R) ("1 L 2 (K). Then by Lemma Ei] 



Az£(|i?/ +Q + |S/ +Q ) <C, 



j 

for any a G [0, 1] and some constant C which is independent of Ax. We also recall 
that for any smooth function / we have that 

if, c j+1/2 D + fj ^'7; (D + R 3 f = cif s (R 2 - Sf) , (A.l) 

where fj = f(Rj), f'j = f'(Rj) and /" = f"(rj) for some rj between Rj and Rj+i- 
We now let a be a constant in [0, 1) and define / to be a C°° function such that 

f(K) -- 



f(K) = I f'(a) 
Jo 




l+a 

Note that f"(K) is bounded. Let x{ x ) be a smooth function such that < x( x ) < 1 
and 

'o, xg[a-l,b+l], 
1, x £ [a, b], 



where a < b are real numbers. Set Xj = x( x j)- We multiply (A.l I by x^Ax, sum 
over j and integrate over [0,T] to end up with 

Ax E &Xi lo - [ Ax E XiCj+ipD+fi dt +£axJ2 Xi^-fj {D+Rjf dt 



( T AxYsXjf&iRt-Sf) dt. 
Jo 
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Figure 2. The scheme ( |4.3| ) (top) and ( |4.4| (bottom), with the 
initial data (4.8) and c given by (|4.7[). 



After a partial summation of the second term on the right, we obtain 
f T 

Ax E fiXi lo + / Ax E fe 2 ^ ^ + c 3-i/* D -Xi] dt 

- [ Ax E ( d + r > f '" A - r E - ^) *• 
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Rearranging this we find that 



JO 



-\R, 

1 + a 



l+a 



(.ft,- - S'j) 

= - Ax hXi \o ~ C I" Ax J2 Xi&i ( R j - S o) dt 
j J ° i 

[ AxJ2 Xfcj ;(//; - S?) ./•; - 2 (//, - 5,-) /,] 1,/, dt 

£ A,-Y.\,^'f';ii)-"f dt. 



By the LP estimates, Lemma |3.4[ all terms on the right-hand side of this are 
bounded by a constant CT, a ,b depending only on T, a, b and on the L 1 and L 2 
norms of Rq and Sq. Furthermore, by the same lemma, 



/ Az^XjSCj 
Jo 



-\Rj 



1 + a 



l{fl 3 >-i} 



Therefore we get the bound 

_,XjCj (Rj-Sj) 

j 



Jo 



\R 



ll+Q 



l + a 



2 ~ 5 i) R J 



dt 



and similarly 

fT 



/ ^ 

Jo 



XjCj 



- (R, - Sj) 



I l+a 



l + a 



dt 



< Cx, a ,b, 

(A.2) 

< Cr,a,b- 

(A.3) 

Adding these two and recalling that = —Rj and |Sj| = — Sj, we get the bound 

^Axj^xfci (ii?,r +Q - i^-r +Q ) m - i^d 
- 1 (i^i 2 - i^i 2 ) ci^-r - i^r) 

The term in the square bracket above can be rewritten as 

o< T ^(i^r +a -i^r +Q )(ii?,i-i^i) 



dt 



<C T , a ,b- (A.4) 



1 

2 

1 1 



\Rj\ -\SjY ){\R j \ a -\s j \ a ) 



l + a 2 
1 



(|^|-|5 J |)(|iZ J | 1+ °-|5 i | 1+Q 



m a \sA a {\Rj\-\s s \) (li^l 1 - -!^! 1 - 



1 



(l-a)(\Rj\-\Sj\) 



2(1 + a) 

x - \s,\ 1+a + \R 3 \ a \Si\ a (\Ri\ l ~ a ~ IS/"" 



■ 2a (\Rj\ - |^|) \Rj\ a l^-r (l^| 1_a ~ ISil 1 "*] 
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' " '\Ri\-\Sj\)*QRj\ a + \Sin 



2(1 + a) 

+ ^ (\Rj\ - \S S \) \Rj\ a \S 3 \ a (l^l 1 - - 



Hence, multiplying (A. 4 1 by 2(1 + a), we arrive at 

i-T 



Jo 



(l-^d^-i-i^ifd^r+^n 



■2a\R j \ a \S j \ a (\R j \ L - a -\S j \ 1 - a 



(A.5) 



dt < C T ,a,, 



Both terms in the sum and integral above are positive, and thus the integrals of 
the sums of the individual terms are also bounded. 
We can use the inequality 

\R\ a + IT > C a \{\R\ - \S\)\ a = C a \R- S\ a 

for some constant C a depending on a, to get the bound 

f AzYVc,- \R, - S 3 \ 2+a dt < C aM (A.6) 
Jo 



Since C\ < c(u) < C2 (cf. ( 1.9 1 ) we have proved the following lemma. 



Lemma A.l. Let a S [0, 1), and assume that Rq and Sq are nonpositive, and in 
L 1 (M) n L 2 (K). Then we have the estimate 

I Ax V c'(u+) \D+u 3 _ 1/2 \ 2+a dt < C a , T , a , b , (A.7) 

where j a Ax £ [a — 1, a) and jbAx € (6, b + 1] . 
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